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ABSTRACT: 

We have shown that photoionization by the continuum of the hidden Seyfert I nucleus in NGC 
1068 can have a significant effect on the ionization state and energetics of this disk's ISM. 
Photoionization models with appropriate power law spectra can produce [NETJ^A. 6538, 6584/Ha. 
line ratios of 1.25 for ionization parameters Q ~10‘ l2 . However our data indicate large regions 
where the [NII]/Ha ratio is 1-3. Since the abundances are known to be solar, there must be 
additional heating sources. Hardening of the incident radiation field by intervening absorption 
should be able to raise Tfe, thereby raising the [NII]/Ha ratio. Heating with moderate efficiency by 
the intense starburst ring should also be a significant factor in raising the temperature of the ISM. 

Our photoionization models with additional heating predict enhanced emission from other 
forbidden lines including [OH] A. 3727 and [SH] A, 6731. 

INTRODUCTION.^ 

NGC 1068 is a prototypical Seyfert II galaxy which has been the subject of considerable study 
due to its close proximity and wealth of interesting activity. There exists in the inner 15".0 a bi- 
polar radio jet, a stellar bar as well as an elliptical ring of intense star formation, with semi-major 
axis parallel to both the bar and radio jets. That NGC 1068 harbors a hidden Seyfert I nucleus has 
been established by the detection of a typical broad-lined spectrum in scattered polarized light 
[Antonucci and Miller, 1985] as well as by the detection of broad Fe II nuclear line emission 
[Snijders, Netzer and Boksenberg, 1986]. Theoretical work [Krolik and Begelman, 1986] 
suggests that the Seyfert I nucleus is hidden from view by an optically and geometrically thick 
toms, which is the source of an electron wind boiled off of its inner surface. Detailed kinematical 
analysis of our Fabry-Perot data indicate that this toms has an opening angle of ~85°, and 
symmetry axis inclined 35° out of the plane of the disk. This implies that the Seyfert I nucleus 
beams substantially into the disk. NGC 1068 has been found to contain considerable quantities of 
spatially extended high excitation gas, particularly toward the northeast quadrant. Evans and 
Dopita, studying HII regions in NGC 1068 [Evans and Dopita, 1986], have determined that their 
abundances are solar and their spectra are superpositions of a normal HII region component as well 
as a high excitation one, characterized by excessive [NeV] X 3346, 3426 and He II X 4686 
emission. 

Our data were obtained using the Hawaii Imaging Fabry-Perot Interferometer (HIFI) [Bland 
and Tully, 1988]. Summations over the separate monochromatic flux images in both Ha and 
[Nil] A. A. 6538, 6584, produce complete Ha and [Nil] flux maps (see Fig. la). A flux ratio map of 
[NII]/Ha (see Fig. lb) is characterized by regions where [NII]/Ha < 0.5, with line widths of 
-100 km/sec, and by kinematically distinct, high excitation regions where 0.5 < [NII]/Ha < 3.0, 
with line widths typically -350 km/sec. 

LARGE SCALE EXCITATION: 

1. Photoionization and heating by a Seyfert I power law continuum 

The 100 eV to 10 keY nuclear spectrum of NGC 1068 can be characterized by a two 
component power law. To determine the shape of the spectrum throughout the disk of NGC 1068, 
we have employed a Monte Carlo calculation to simulate the scattering of radiation out of the 
beaming cone over the energy range of 7.5 eV - 1 MeV. The direct radiation within the beaming 
cone, as well as that scattered into the disk outside of the cone are given in Table 1. Due to the 
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energy independence of low energy Compton scattering the spectral index for the 0-5 keV 
component does not change. At higher energies the energy dependence of both the cross-section 
and the energy transfer to die electron tend to steepen the spectral slope. 


ENERGY RANGE 


TABLE h 
SPECTRAL INDEX 


FLUX AT 13.6 eV 
(ergs/cm 2 sec Hz) 


Within Cone 


0 - 5 keV 

2.0 

2 x 10' 15 (1 kpc/D) 2 

5 keV - 1 MeV 

0.62 

lx 10-18 (lkpc/D) 2 


Outside Cone 


0 - 5 keV 

2.0 

1 x lO' 12 (1 kpc/D) 2 

5 keV - 1 MeV 

0.81 

5 x 10-21 (\ kpc/Dl 2 


We have used the general purpose ionization code MAPPINGS [Evans and Dopita, 1985] to 
investigate the ionizing effects of such an input spectrum. We have modeled the photoionization 
by a geometrically thin spherical shell with volume filling factor 0.1 and density 1 cm' 3 . Since the 
code has a high energy maximum of 5 keV we have chosen a spectral index of 2.0. We find that 
we are able to generate [NII]/Ha ratios of up to 1.25 for ionization parameters Q = 10" 12 where Q 
is defined by 


Q= 


_ 1 _ 

NtotC 


r 


dv 

hv 


However these models can not simultaneously generate the high values of [NeV] XK3346, 3426 
and Hell X 4686 emission seen in the Evans and Dopita data. These values require ionization 
parameters an order of magnitude greater or more, indicating that these emissions may be spatially 
distinct. 


We have calculated the Compton heating caused by photons of energies greater than 5 keV. 
For typical cooling rates for ionized gas of 10' 24 N e 2 ergs/cm 3 sec, we find that the additional 
Compton heating in the cone will only be significant for ionization fractions less than 10' 1 and/or at 
very low densities, whereas the Compton heating off of the cone should be negligible under all 
realistic physical conditions. 

2. Heating by the Starburst Ring: 

It is difficult to determine the direct heating effect of the starburst ring on the extended disk of 
NGC 1068 due to the undetermined geometric and starburst parameters. However an estimate may 
be made by scaling up the starburst of M82 which is better understood. Using a scaling factor 
equal to the ratios of the IR luminosities of these two starburst systems, we can determine the SN 
rate in NGC 1068 to be = 0.75/yr, which can generate = 2 x 10 43 ergs/sec. Assuming a 3 kpc 
radius and 1 lq)c thickness for the hot disk of NGC 1068 with a filling factor of 0.1 and mean 
filament density of 1 cm -3 , we need ~ 1% efficiency in converting SN energy into heat in the disk, 
to have a heating rate of 10' 24 ergs/cm 3 sec, comparable to that of canonical heating/cooling rates. 


144 




Fig. la. Total Ha Flux Map Fig. lb: [ND]/Ha Ratio Map 


REFERENCES: 

Antonucci, R. R. J., and Miller, J. S., Ap. J., 297, 621, 1985. 

Bland, J., and Tully, R. B., Astron. J., August, 1988. 

Evans, I. N., and Dopita, M. A., Ap. J. Suppl.., 58, 125, 1985. 

Evans, I. N., and lopita, M. A., Ap. J. Lett., 310, L15, 1986. 

Krolik, J. H., and Begelman, M. C., Ap. J. Lett., 308, L55, 1986. 

Snijders, M. A. J., Netzer, H., and Boksenberg, A., MJSI.RA.S., 222, 549, 1986. 


ORIGINAL, PAGE IS 
OF POOR QUALITY 


145 





